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We present a Nuclear Magnetic Resonance experiment aimed to scan the evolution of a spin sys¬ 
tem, in liquid crystals, during the transient to a quasi-equilibrium state. New evidence is presented 
in favor of irreversible decoherence as the mechanism which leads an initial out-of-equilibrium state 
to quasi-equilibrium. The experiment combines the Jeener-Broekaert sequence with reversal of the 
dipolar evolution, and decoding of multiple quantum coherences, to allow visualizing the evolution 
of the spectra of the different coherences during the formation of the quasi-equilibrium states. We 
vary the reversion strategies and the preparation of initial states, and observe that the amplitude 
of the spectra attenuate with the reversion time, and notably, that the decay is frequency selective 
in all the samples. We interpret this effect as an evidence of “eigen-selection”, a signature of the 
occurrence of irreversible decoherence, which indicates that the spin system in liquid crystal NMR 
experiments conforms an actual open quantum system, where quasi-equilibrium can be rigorously 
described as a stage of the dynamics which develops in a time-scale far earlier than thermal equi¬ 
librium. This explanation is supported by a derivation of the observed signals and spectra in terms 
of the coherences of the spin-system of a single molecule, within the theoretical framework of open 
quantum systems. 

I. INTRODUCTION 

The quantum dynamics that drives a quantum system 
of interacting particles from an initial coherent state into 
a thermodynamic state are of high concern in physics. 

The interest relies on fundamental questions as well as 
on practical applications. Nuclear spins are privileged 
systems whose quantum dynamics develops in a wide 
range of time scales that can be accurately probed with 
Nuclear Magnetic Resonance (NMR). The possibility of 
preparing different initial quantum states in diverse spin 
systems, turn them into excellent test-beds for studying 
the interesting quantum dynamics during their transit 
from an initial coherent state to the final thermal equi¬ 
librium. Particularly, our interest in this work is focused 
on the intermediate stages, where the irreversible trans¬ 
formation of the initial coherent state into the so called 
“quasi-equilibrium” (QE) states takes place. Starting 
from a coherent initial state, QE is attained over a tran¬ 
sient period, longer than the time scale characterized by 
the typical spin interactions. The spin observables as¬ 
sociated with this state evolve subsequently only due to 
spin-lattice relaxation, within a much longer time scale. 

The mechanism governing the build up of the QE states 
is still not fully known. It is the aim of this work to 
provide experimental evidence on its origin. 

We consider the spin system formed by the proton 
spins in liquid crystals (LC). In this kind of sample each 
molecule constitutes a small cluster of interacting parti¬ 
cles, coupled to a highly correlated non-spin environment. 

In fact, due to the rapid individual molecular motions, 
the intermolecular spin interactions average out, still, a 
strong intramolecular residual dipolar energy remains be¬ 
cause of the high degree of long range molecular orien¬ 


tation typical of the mesophase m- The NMR spec¬ 
troscopic properties of LC are similar to those of solids, 
that is, the on-resonance free induction decay (FID) sig¬ 
nal evolves over a short time scale determined by the local 
dipolar couplings in the molecule, and the spin-lattice re¬ 
laxation time scale is much longer than that of the FID. 
Likewise, it has been shown that in spite of the small 
number of relevant spin degrees of freedom, the LC spin 
system can also develop quasi-equilibrium states within a 
time scale similar to solids nig. Indeed, the experimen¬ 
tal results clearly showed that it is possible to selectively 
prepare several QE states in LC by adequately setting 
the preparation pulse sequences liE]. 

Owing to the small size of the spin system, the at- 
tainance of QE in LC cannot be explained through ar¬ 
guments borrowed from statistical mechanics, based on 
the occurrence of a large number of interacting spins, like 
spin diffusion. In previous work we anticipated that the 
development of QE should be related with irreversible 
adiabatic decoherence ElE!, that is, with the environ¬ 
ment induced decoherence without dissipation. Such a 
process erodes the off-diagonal elements of the spin den¬ 
sity matrix in the “preferred” basis m, driving the sys¬ 
tem to a state with a diagonal-in-block form. 

We studied the irreversible attenuation of the NMR 
signals in refocusing experiments (MREV [3] and Magic 
Echo m), and observed that the irreversible decoherence 
and the development of QE take place within the same 
characteristic time scale. That is, over an intermediate 
time scale, considerably earlier than the one dictated by 
thermal processes. Further theoretical work liiaiii] in¬ 
dicated that a small cluster of interacting spins, viewed 
as an open quantum system coupled to a correlated en¬ 
vironment, can indeed undergo irreversible adiabatic de- 
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coherence, just as the experimentally observed in LC. 
A very important feature predicted by this theoretical 
proposal is that decoherence is controlled by the eigen- 
frequencies of the spin part of the interaction Hamilto¬ 
nian rather than those of the system Hamiltonian. This 
phenomenon implies a kind of selection (induced by the 
coupling with the environment) since each element of the 
density operator decays with a characteristic time, scaled 
by the corresponding eigenvalue differences, and was ac¬ 
cordingly called “eigen-selection” [12] . The occurrence of 
this effect was experimentally observed on the single-spin 
single-quantum spectra of some nematic LC [S] . The ex¬ 
periment showed that the higher frequency components 
of the spectra of the reverted signals attenuate with a 
higher rate than the lower frequency components, and 
this behaviour cannot be attributed to experimental non¬ 
idealities (non-reverted terms, pulse setting errors, etc). 

In this work, we present an NMR experiment designed 
to display the evolution of coherences during the buildup 
of a QE state. The technique is applied on several liquid 
crystals. The experiment aims to visualizing the spectra 
of the different coherences during the formation of the QE 
and combines the Jeener-Broekaert (JB) sequence [13] 
with reversion of the dynamics and encoding of all the 
coherence orders. This method allows to show the effect 
of eigen-selection during decoherence, which produces the 
characteristic frequency compression on the spectra. We 
are particularly interested on exploring the transient to 
QE states other than the traditional dipolar order (two- 
spin order), featured by a higher multi-spin character. 
We expect to detect the traces of decoherence in the non¬ 
equilibrium MQC spectra. An additional convenience of 
the small-sized spin systems of LC is the possibility of 
performing the exact numerical calculation of the den¬ 
sity operator that would be expected in a closed-system 
evolution in a small spin cluster. This allows comparison 
between the calculated and experimental NMR signals. 
In this way, we bring new evidence which suggests that 
decoherence in spin systems is produced by the quantum 
coupling with the environment. 

In Section |ll| we explain general aspects of the experi¬ 
ment which tracks the evolution of coherences and spec¬ 
tra during the buildup of quasi-equilibrium states and 
Section [ha contains an analytical description of the sig¬ 
nals. The particular case of the output spectra in ne¬ 
matic liquid crystals treated as open quantum systems is 
described in Section |HB Section shows the experi¬ 
mental results, and in Section |TV] we expose the conclu¬ 
sions of the work. Appendix A is added to show the effect 
of non-idealities on the outcome of the experiment on a 
hypothetically closed system. 


II. DESCRIPTION OF THE PULSE SEQUENCE 


The experiment was tailored to study the way in which 
decoherence affects the frequency content of the spec¬ 
tra of different coherences through eigen-selection [El IS] 


during the transit from a multi-spin initial state to the 
quasi-equilibrium. The sequence combines the Jeener 
Broekaert (JB) experiment with encoding of multiple 
quantum coherences and refocalization of the spin dy¬ 
namics under the secular dipolar Hamiltonian. 

In the standard JB sequenc^(7r/2)2, —tp —(Tr/J)^ —t — 
(7r/4)p+,p] (top row in FigurejHa), with r = 0), the first 
(7r/2)a; pulse produces an initial state of one-spin single 
quantum coherence which evolves freely under the secu¬ 
lar dipolar Hamiltonian during tp, developing more com¬ 
plex multispin, single-quantum coherences. The second 
pulse (7r/4)p rotates the state into a multi-spin multiple- 
quantum state (zero, single, double or higher order). The 
state so prepared is allowed to evolve freely during a 
waiting time t. If t is “long enough” the spin system 
reaches a QE state consistent with a diagonal density 
operator (in the basis of the dipolar Hamiltonian); let 
us call tQE to the time scale needed for the particular 
spin system to reach quasi-equilibrium. A last read pulse 
(7r/4)p projects part of the QE state into observable sin¬ 
gle quantum coherence. The QE state which can be ob¬ 
served after the read pulse depends on the preparation 
time, in fact, it was shown that at least three different 
QE states can be observed in nematic liquid crystals [7]. 

The experiment proposed in this work, depicted in Fig¬ 
ure Sa) , aims to make visible the effect of decoherence 
during the development of quasi-equilibrium, we thus (i) 
prepare multiple quantum coherences, {ii) display their 
spectra, and (Hi) observe the spectra during their evo¬ 
lution towards quasi-equilibrium. We prepare the initial 
state with the first phase shifted pair of pulses of the JB 
sequence and use the phase ip of the read pulse to encode 
the different coherence orders M- The experiment is 
repeated while incrementing the waiting time t in steps, 
so that a Fourier Transform on t and ip gives the spectra 
of each coherence order separately |I5j . The important 
information of the effect of decoherence on the spectra is 
probed by adding a block, marked as D in FigurejHa), 
at the end of which the free evolution under during 
time T has been reverted. The effect of block D is thus to 
“filter” the reversible evolution. The consequences of the 
irreversible mechanisms on the signal are then observed 
by incrementing the duration of block U, as shown in the 
rows in Figure [^a). 

We use two different schemes for the reversed dynamics 
to ensure that the results are independent of the partic¬ 
ular experimental strategy. Some compounds were stud¬ 
ied with the MREV8 sequence p iTHHIE] . schematized in 
Figure [^b) and other with the ‘magic-sandwich’ (MS) 
sequence [TTJ[TnHE3] (see Figure [^c)). Each sequence of¬ 
fers particular advantages: MREV can be both used with 
a continuous variation of the time ti or by concatenating 
blocks of this sequence [9] , so allowing a thorough scan of 
the decoherence time; MS involves less (although longer) 
pulses, and is consequently less vulnerable to errors in 
the pulse settings. 

In summary, Figure [^a) illustrates the whole exper¬ 
iment, where a set of sequences at different refocusing 
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FIG. 1. (a) The pulse sequence used to trace the formation of 
quasi-equililibrim starts with the phase-shifted pulse pair of 
the JB sequence, where tp serves to select the QE state that 
will develop at the end of the sequence. The waiting time 
t is incremented in steps as well as the phase ip of the read 
pulse. Fourier Transformation on t and 93 gives the spectra of 
each coherence order separately. Within block D of duration 
r, the system undergoes a reversion dynamics while evolving 
towards the quasi-equilibrium. The reversion strategies used 
in the experiment consist of (b): the MREV 8 sequence, where 
r = 12 ri; or (c): the ‘magic-sandwich’ sequence, where r = 

1.5 tm 


times r are vertically arranged in order to imitate the 
bidimensional disposition chosen to show the experimen¬ 
tal results (see Section III). Each row represents the set 
of Ni X experiments needed to calculate the spectra 
of all the encoded coherences (where Nt is the number of 
steps of t and the number of phase steps in the read 
pulse). The next sections contain a derivation of the out¬ 
put signals of this experiment and the interpretation of 
the corresponding spectra in the framework of the theory 
of open quantum systems. 


A. Description of the signals and spectra 


When the observed spin system in an NMR experiment 
is considered as an open quantum system in contact with 


an external bath, the total Hamiltonian can be written 
as 


n = ns + nL + nsL, (i) 

where the system of interest is described, in isolation, by 
Us] 'Hl the Hamiltonian of the external reservoir and 
represents the system-bath interaction which has 
both spin and non-spin variables (e.g. interaction be¬ 
tween spins and mechanical variables, or a boson bath, 
etc.). The hat symbol over operators stands for the “ro¬ 
tating frame” representation; one which rotates around 
the quantization axis z at the Larmor frequency, and 
where the spin Hamiltonian "Hs does not include the Zee- 
man term Hz = — wo Iz- 

Since we are interested in describing processes which 
take place in a timescale earlier than relaxation and ther- 
malization effects, the Hamiltonians satisfy the adiabatic 
condition [51 IMll^ 


Hs,Hsl 


= 0 . 


( 2 ) 


Assuming adiabaticity is equivalent to excluding en¬ 
ergy exchange between the spins and the bath and allows 
factorizing the total evolution operator as 

Ur(t) = ‘e-* t ^ iji^t)ijsLit), (3) 

where the evolution operator 

U(t) = e-*^«‘ (4) 

acts on the spin degrees of freedom only and describes the 
dynamics of a closed system. The “openness” enters in 
the evolution under the non-commuting terms Hsl+Hl 

fjsL(t) = (5) 

See Section H 2 from [5] or Section H A from [5] for 
details about the dynamics under these operators. 

The initial equilibrium state of the complete system 
can be written as the tensor product P(^eq) = Ps{eq)PL(eq)^ 
where Ps{eq) oc Iz represents the spin system in equilib¬ 
rium with an external field Bq in the high temperature 
regime, and PL{eq) is the density operator of the non¬ 
spin degrees of freedom (which is not affected by the RF 
pulses nor by the evolution under only-spin operators). 

The first {tt/2)x pulse leaves the spin state as 

fc(0) = Rx(7r/2) I, Rx(-^/2) = ly, (6) 

where the pulse operator Rq( 0) = represents a ro¬ 
tation through an angle 9 around the direction a. 

After the first pulse, the state evolves under the total 
evolution operator of Eq. ([^. However, it is worth to 
notice that both the preparation time tp and the acquisi¬ 
tion time t' , used in the experiment outlined in Figure [l] 
are much shorter (few tens of microseconds) than the 
timescale scanned by the sum t + t (tens to hundreds 
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of microseconds). Since the dynamics driven by Usi, is 
slower than that driven by U, we can consider that 

Usz,(t) ^ 1 during short times t, (7) 

which justifies assuming a purely unitary dynamics 
during tp and t'. Then, under the assumption Q 
the state at the end of the preparation period is 
p(tp) = U(tp)ps(0)Ul(tp) On the contrary, we 

expect the effect of decoherence due to the inevitable quan¬ 
tum coupling with the environment, to manifest during 

T + t. 

The second pulse (7r/4)j, rotates the spin state devel¬ 
oped at tp and transforms it in a spin density operator 
with terms of zero, single, double or higher order coher¬ 
ence, and the state immediatly after preparation is 

p(t+) = Ry(7r/4)ps(tp)Ry(-7r/4) (8) 


I(x,y)-|-(p — R'z( T’) I(x,y) R-z(<p)! (12b) 

and after some algebra, the signal can be written as 


Sa = tr{laij{t')'Ry{TT/4:)'Rff(p)a{t,T,tp) 
X R,(-^)Ry(-V4)Ul(t')}. 


(13) 


In order to disclose the coherence encoding implicit in 
Eq. (131, it is useful to write the reduced density opera¬ 
tor in terms of a set of irreducible tensors which 

serve as an operator basis for the Liouville space where 
a belongs. Index A is the tensor rank, v the component 
number, and A distinguishes between different species of 
tensors with the same rank m- The reduced density 
operator at any time is formally 


r{t,T,tp) = 


,A 


(14) 


A \v 


After the second pulse, the state is driven through a 
forward-backwards evolution during r, which has the ef¬ 
fect of undoing the evolution under the secular dipolar 
Hamiltonian, but, in principle, cannot undo the evolu¬ 
tion under fisL + ‘Hl- Let us represent this evolution 
by the operator U_r(t) [ 5^. Then, the density opera¬ 
tor after the reversion period (r) and the free-evolution 
(during t) is 

p{t,T,tp) = Ut( 1) Ufl(r)ps(t+)pi(e,)U]^(T) Vl^it) . 

(9) 

Since this state is no longer separable and we are con¬ 
cerned in calculating the signals produced by the spin 
observables, we use the spin reduced density matrix de¬ 
fined as 


a{t,T,tp) = tr^ {p{t,T,tp)}, 


( 10 ) 


where the partial trace tr^ {•} runs over the environment 
variables. It is worth to note that a as well as ps are both 
density operators which apply only over the spin Hilbert 
space, but the definition (10) implies that the dynamics 


that affects to a can be non-unitary. 

Finally, a read pulse is used to rotate the tensor com¬ 
ponents with arbitrary v, into one with v = ±1. The read 
pulse phase tt j2 -\- ip (represented as y -|- in the pulse 
operator) is varied systematically (and consequently the 
receiver phase) in order to encode the coherence order. 
Therefore, the NMR signal S, as a function of the acqui¬ 
sition time t' and the other variable parameters of the 
experiment {p,t,T and tp), is 


Sa{t',g},t,T,tp) = tr{la U(t') R-y-Hv(^/4) 

X CT(t,r,tp)Ry+,^(-7r/4) V\t')}, 


( 11 ) 


where a = x or y. In Eq. (11) we also used the approxi¬ 
mation ([^ because we only need the signal at the earliest 
acquisition times t' . Using that 

Ry+<p(7r/4) = Rz(-(/?) Ry(7’'/4) Rz(‘/7), (12a) 


where 

CA(l>Ulp) = tr \T{^a{t,T,tp)'^ 
is the projection of a{t) on which satisfies 

and 6 is the Kronecker delta. 


(15) 


Finally, using expansion (14) in Eq. (13), and the fact 
that 


Rz(</p)TARz(-7’)=e*'^^Ti, (16) 

the NMR signals measured under the sequences of Fig¬ 
ure can be expressed as 

Sa{t',ip,t,T,tp) ='^e^'"'^'^gxl°'{t')f.t{t,T,tp), 

V A A 

(17) 

where 

= lr{l„U(T)Ry(7r/4)TARy(-7r/4)Ut(t')}. 

(18) 

In our experiment, instead of we used g^)^{t'^), 

the average value of the signals within an interval A cen¬ 
tered at the earliest measuring time t' = t'^ where the 
signal has its maximum absolute value in some experi¬ 
mental outcome. That is 

5t“(C)^x/ (19) 

^ Jty-Al2 

The spectra of coherences emerge after calculating the 
Fourier transform of Eq. 0 over the variables t and tp, 
that is 

•A’l* (im! Ip) “ 

n)'^g^l°‘{t'm)J^t {CA} (w,T,tp), (29) 

ly A A 
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where Ta {/(o^)} (/3) represents the Fourier transform of 
function / over the variable a, and /3 is the conjugate 
variable of a. In Eq. (20), /r and w are the conjugate 
variables of and t, respectively. 

Equation (20) stands for the coherence spectra ob¬ 
tained at a given combination of preparation time tp and 
evolution under reverse dynamics in r. The novelty it 
introduces comes in the non-unitary dynamics along the 
intermediate time scale r, which enters in the observed 
spectra through the Fourier transform of the coefficient 
^{t,T,tp). If the observed system were actually closed, 
its NMR signal would not depend on t because the dy¬ 
namics would be completely reversible. However such 
a reversible behaviour could be a fictitious expectation 
since actual physical systems are not closed, indeed, the 
subtle quantum correlation between the observed system 
and its environment (even without energy exchange) may 
cause decoherence [an]. Thus, any observed r depen¬ 
dence might be interpreted either as an indication of an 
irreversible (non-unitary) behaviour, or as the result of 
an incorrect experimental setting. 

In the next section, we write (20) for the case of ne¬ 
matic liquid crystals (NEC), where the eigen-selection 
effect on the spectra of the coherences under reversion 
dynamics, becomes evident mm- 


B. The case of nematic liquid crystals 


This section is dedicated to obtain an analytical ex¬ 
pression of the output signals of the experiment of Sec- 
tion|^on a NEC sample, considered as an open quantum 
system. We first analyze the reasonable assumption that 
in a state prepared at tp there is no particular correlation 
among the spin variables of different EC molecules. How¬ 


ever, the subsequent evolution under UT(t) and Ufl(r) 
may mix the Hilbert spaces of different molecules and 
turn the state into a correlated one. Then, in order 
to facilitate the analysis of the experiment in terms of 
the charcteristic features of the nematic ordering, in Sec¬ 
tion HB 2 we write the state at the end of the pulse 
sequence described above in terms of the state of a single 
molecule. 


1. Initial non-correlated condition 

The main contribution to the total spin Hamiltonian in 
NEC comes from the dipole-dipole interaction. The fast 
molecular diffusion occurring in the nematic state, av¬ 
erages out the inter-molecular dipolar interaction, leav¬ 
ing only the intra-molecular part [2] . This characteristic 
of orientationally ordered phases allows writing the spin 
Hamiltonian Hs of a NEC in a separable form mm , in 
terms of the single molecule Hamiltonians 

= ( 21 ) 


with - - where the su¬ 

perscripts between parentheses indicate different Hilbert 
spaces: (s) for the whole spin sample, (s^) for the spins 
that belong to the i-th molecule and (e) for the envi¬ 
ronment, non-spin variables. is the secular dipolar 

Hamiltonian of a single molecule, in the high-held ap¬ 
proximation 


E/l®*) — ^ c 

ric — 


(i) 






( 22 ) 


The symbol (22) represents a sum that 


runs over the sites of the resonant nuclei within the i- 
th molecule, and is the secular irreducible spher¬ 

ical tensor of second rank which applies over the Hilbert 
space of the i-th molecule, 


2Iz)/1zA: T I—jl-t-fc) j 

(23) 

where lak is the a-component of the angular moment of 
the fc-th spin and indicates that the operators apply 
exclusively over the Hilbert space of the z-th molecule. 
The dipolar frequency dehned in Hz units is 

ujoirkj) = 3^^^ [1 - 3cos2(/3jfe)] , (24) 

with j3jk the polar angle of vector joining spins j and 
k, with respect to the z-axis of the system fixed to the 
molecule and is the mean value of the molecular order 
parameter |S], defined as 


r^jk{si) 

^20 


76 





(25) 


with 9 the polar angle between the molecular z-axis and 
the static strong magnetic field Bq. 

A consequence of the separability of Hamiltonian (21) 
is that the inter-molecular correlation is not involved in 
the Eiouvillian evolution, within the adiabatic time scale. 
Consequently, the evolution operator (i can be factor¬ 
ized as 


U(t) (26) 

with - • •(g)0('’‘)(g) - • - ^where it 

is clear that U (t) acts on the spin space of each molecule 
without mixing their Hilbert spaces. 

Eet us now consider that the experiment starts from 
thermal equilibrium. Therefore, the spin part of the den¬ 
sity operator is 


-(s) _ 6' 

Psie,) = - 


./3t 


l(s) 

K 




PsLy ( 27 ) 




where 0 ••• 0 0 ••• 0 and 

A/si = tvsi is the trace of the identity operator 
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over the Hilbert space of spins belonging to only one 
molecule, N is the number of molecules, and we define 
the equilibrium density operator of the single molecules 





(28) 


with fix = the Boltzmann constant and T the 

absolute temperature. The complete thermal equilibrium 
density matrix, that involves spin and environment vari¬ 
ables, is 


'^{s) „ (e) 

P{eq) — Ps(eq) ® PL(eq)- 


(29) 


In this way, using Eq. (26) and the initial thermal state 

immediately 


(27), the spin density operator from Eq. 


after the preparation period, for a NEC can be written 
as 




(30) 


Si 


where p_^®^(f+) = 0 • • • (g) p(®‘)(f+) (g) • • • 0 and 




■^Si 


Ry(V4)U(tp)I, 


xUt(gRy(-7r/4) 


(S.) 


(31) 


It should be mentioned that we leave out the identity 
operator in Eq. (30) because it does not contribute to 
the NMR signals. 

Therefore, the complete density matrix for the state is 
++^ — ^ (32) 


m)=n’it;)^P'Lieqy 


In obtaining Eq. (32) (and thus in Eq. (30) as well), we 
used the approximation 0 and Eq. ([^. 

Finally, Eq. (32) is the non-correlated initial state 
whose dynamics under decoherence we calculate in the 
next section. 


2. Evolution of the spectra in an open quantum spin system 


We now calculate the dynamics under decoherence 
in the same way as in Eq. starting from an non- 
correlated initial state like (32). In this work we will use 


the expressions for the evolution operators obtained from 
reference [5] , which describe the non-unitary dynamics of 
a NEC sample [2H1. 

In order to visualize the coherence content of the 
observable signal as in Eq. dTi 


and (20), we follow the 
reasoning which goes from Eq. ([^, defines the reduced 
density operator (10) and spans it on the tensor basis 


of Eq. (14). It is worth to have in mind that within 


the time scale where coherences may grow and atten¬ 
uate, the evolution operator (|^ introduces correlations 


between the spin Hilbert spaces of different molecules by 
means of the non-spin or environment variables. There¬ 
fore, the reduced density operator (which by definition 
applies exclusively over the spin space), cannot in princi¬ 
ple be written as a linear superposition of tensors acting 
independently over each molecule. However, we can de¬ 
fine a tensor basis for the complete spin Hilbert space as 
the tensor product of the elements of the single-molecule 
bases, In this way, we can write the ex¬ 

pansion of the reduced density operator at any time t 
as 


aW(t,r,tp)= ••• II ••• H 

AiAii/i A-iXfi^i 


(33) 


with 




-ifAifsi 
■ Xii'i 


(s) 


{t,T, tp)| 


CsX's' 

(34) 

where the trace runs over the total spin space. In 
Eq. (34) the indices in curly brackets indicate depen¬ 
dence on the set of eigenvalues of each molecule of the 
whole sample, denoted as {Ai} = {Ai,..., A^,..., Ajv} 
and {Xi^i} = {XiUi,..., XiVi,..., Xjqvjq}. Besides, we 
have the following orthonormality condition 

trs. ^ 6^, 

(35) 

with and = 0. 

In order to calculate the trace in Eq. (34), we use 
a product basis of the common eigenstates of the spin 
Hamiltonian, 'Hs^ and the spin part of the spin-lattice 
Hamiltonian of each molecule (remember that the 

commutation relationship 0 implies the existence of a 
common eigenbasis), that is 

{|Cs) = «)i|Csi)}, 

where Q are the eigenvalues, Si accounts for the de¬ 
generation of each eigenvalue and the eigenstates satisfy 

= S'zzCil CiSi)- 

The elements of the reduced density matrix of Eq. (34) 
are 

(Cs|CT^'^Ht,r,tp)|C's') = (Cs|tre|p(®)(t,T,tp)| |C's') 

i ^ Si 


^ If ^{C.C'}(^>^)- 




(36) 
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where (t, r) is the decoherence function associ¬ 

ated with the non-unitary evolution due to the spin- 
environment coupling. This function is defined from the 
general expression 

(37) 

where the indices C and C represent the set of eigenval¬ 
ues of all the molecules, that is ^ = {Ci,..., Ci) • ■ • > Cn} 
and ..., ..., The evolution operator 

f, r) acts on the environment degrees of freedom 
only, and is defined from the evolution operators Ur(t) 
and Ujj(t) (which apply over the spin and the environ¬ 
ment Hilbert space) as follows 


function and also to exhibit its relation with the observed 
signals and spectra, we first assume the existence of dif¬ 
ferent time scales for different mechanism of decoherence. 
As can be seen from Section II B of reference [3] , the de¬ 
coherence factors that affect to the f-th molecule and that 
depend on the eigenvalues (j (with j ^ i) of the rest of 
the molecules have slower decay than that produced by 
the reversible AQD (which depends on the eigenvalue Q 
only), thus they can be neglected in front of the reversible 
decoherence factor. Therefore, the decoherence function 
depend on the spin variables of each molecule only, that 
is, 

= G{(;,,^/}(t,T). (40) 


(C, t, r) = ( Cs I e* * I Cs ) ( Cs I UT(t)Ui^(r) | Cs ) ■ 

(38) 


In this way, the decoherence function in Eq. (36) is ob¬ 
tained by applying the condition Q = Vj/J ^ * on 
Eq. (37). Therefore, depends on all the values 

Qi,. ■ ,C,M and on Eurther properties of the 

decoherence function are explained in reference [3]. The 
product symbol used in Eq. (36) stands for the product 


of the iV — 1 Kronecker deltas of all the molecules except 
the one with j = i, that is 






Cjv^iY’fw^JV ■ 




Using Eq. (36) in Eq. (34), and after of some manipu¬ 


lation with the ordering of indices, we obtain 

i CiS^Xi^i 


w-AT-l 


E 




ITc, 


Aj (sj ) I 


C-5 


J J 


,) 


j¥=i 


(39) 


where the sum ^ runs over the values of Ck^k 

different from (C^Sj. Besides, the sum affects the factors 
in the product and the decoherence function as well [33] . 


The coefficient (39) introduces the effects of decoher¬ 


ence in the dynamics of the density matrix (33), through 


the decay function G{^_^'}(t, r). It is important to re¬ 
mark that this decoherence function emerges in the spin 
density operator as a consequence of recognizing the cou¬ 
pling with the environment and involving it into the spin 
dynamics. 

It is worth to note that, in principle, the decoherence 
function introduces quantum correlations between the 
molecules of the NEC, even when there is no direct inter¬ 
action between spins belonging to different molecules. In 
order to get insight on the significance of the decoherence 


We also make the following plausible assumptions 
(called HypoG-I and HypoG-II in Ref. i)- One is 
the statistical independence between the eigenvalues of 
the spin-bath interaction of each molecule and the eigen¬ 
values of commutators of such Hamiltonians, which al¬ 
lows to separate the decoherence function in two factors 

G{c.,c'}(i,T) = G[^^,^,}(t)Gf^^,^,j(T). (41) 


The first factor represents the strictly adiabatic quantum 
decoherence (AQD) |S], which is responsible for the line- 
shape of the coherence spectra and is the main decoher¬ 
ence process in a free evolution dynamics. It takes place 
in a time scale shorter than the rest of the decoherence 
processes and its dynamics is in principle reversible. The 
second factor accounts for the essentially adiabatic quan¬ 
tum decoherence |8] , it is of irreversible character and has 
a slower dynamics than the reversible AQD. The second 
hypothesis is to consider an homogeneous environment 
for each molecule, i.e. absence of border effects, where 
each molecule sees the same neighborhood, which im¬ 
plies that the function is the same for each molecule, this 
i® = G{c^,q}{t,T) Vi, j. 


In order to simplify the coefficient (39), we see that 
assumption (40) implies that the decoherence function 
in the sum only depends on the spin variables Q and Q' 
of the z-th molecule, and thus can be extracted out of the 
sum s /s which excludes the i-th molecule. We 

thus denote the resulting sum as 


A 


f^J ^3 } 


E 

Cfc 




o¥=i 




(42) 


N-l 


Because of the orthogonal property (35) and the fact 
that the identity operator forms part of the tensor basis. 

















the rest of the basis operators have null trace. Therefore, 


A 




= < (If Vj Aj = 0, Xj = 0 and vj = 0) 


0 (otherwise) 


n 


r0,0,0 

"A,,A, 


,/v^ 


N-1 
Si ’ 


(43) 


where = Sa^,o Sx,,o 

Finally, the commented considerations over the decoher¬ 
ence function together with (43) in Eq. (39), allow writing 






(44) 


where 




X ( C.S. |P C's'z )( C's' )> 


(45) 


which only depends on the coefficients A^, Xi and Vi of the 
j-th molecule, with Tip) = Vi, j. 

Let us now go back to the spin reduced density matrix 
of Eq. (33) and replacing the coefficient (44) we get a 


factorized form 


?(« 


{t,T,tp) = Y 


l(si) 


Si 


tp){ 


]_(««) 

®A5r' 

(46) 

where each factor corresponds to a 

single molecule, and also satisfies a^^'\t,T,tp) = 

a^^A(t,T, tp) MiA- 

It is worth to note that the product of Kronecker delta 
functions that arise on Eq. (44), make the expansion 


(33) to adopt a non-correlated form in the spin vari¬ 


ables, that is, the sum over the indices Xi , Vi and 
of the i-th molecule involve tensor operators which ap¬ 
ply on the Hilbert space of such molecule (with identi¬ 
ties on the indices of the other molecules). We could 
then describe the time dependence of the reduced den¬ 
sity operator of Eq. ( [44| in terms of single-molecule deco¬ 
herence functions which however enclose quantum inter- 
molecular correlations. That is, there is no direct spin 
interaction between different molecules but they become 
correlated because of the more subtle quantum interplay 
between Hs and Hsl 


As a final step, we use the density operator (46) in 


Eq. (13) to calculate the expression for the signal 0 

= NYA''^^ Y 9Yy')^Yi(i^r,tp), (47) 


Ai Ai 


where, similarly to Eq. (|T8|), 

gYyi^tr^i 


I„U(t')Ry(7r/4)T(^),^ 

X Ry(-^/4)Ut(t') 


(48) 


Equation ( |47[ ) shows that the signal generated by the 
total sample of a NCL is N times the signal generated 
by one molecule (here labelled as molecule 1), with N 
the total number of molecule in the sample. 

Finally, the evolution of the spectra for a NLC, mea¬ 
sured by means of the experiment shown in Figure is 

■7)^,4 {'5} (^m5 Ml T tp) 

Y Va( 4)?(0 At {Ca(^^ } (w, r, tp), 


Ai Ai 


(49) 


with gxli^it'm) defined as in Eq. (19). Besides, using 
Eq. (45) in Eq. (49), the Eourier transform in t can be 


written as 
At 


{CA^i}(^.Tip)= ^ Gfci.C(}(^) 

271 /a;-(C(-Ci)^4 


""iCi-CirV C(-Ci 

X (CiSi 

where we use the result 

271 /'uj-{Ci-Ci)Sz 


(50) 


IC(-Ci 


p 


a-Cl 


(51) 


which was calculated in Section II B of reference [9]. In 
( [50| ) the function p is the called orientational molecular 
distribution function (OMDF) [5], which is the Fourier 
transform of the reversible AQD function. 


To conclude this section, we can observe from (49), 


(50) and (51) the following properties: The reversible 


AQD gives rise to the line-shape of every coherence, as 
comes out from Eqs.(|4^ and (1^. Such spectra are ob¬ 


tained as a superposition of copies of the OMDF shifted 
to the frequencies (Ci — Ci)Szz and scaled by the factor 
I Cl — Cl I j as can be seen in ( [m| ). The consequence of this 
scaling is the eigen-selection effect in the time domain 
of t, due to which the larger the value of |Ci ~ CiL the 
faster the decay. This is reflected in the spectra making 
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that for higher frequencies (Ci — Ci)'S'zz the correspond¬ 
ing copy of the OMDF is less intense and wider, which 
produces a resulting effect of narrowing or compression 
of the complete spectrum. 

On the other hand, the irreversible decoher¬ 
ence function j (r) has a decay time which 

depends on the difference of the eigenvalues 
Cl — Cl as well (as an example from reference [S], 

^,{T)=e ^ J ), which intro¬ 

duces an alteration of the line-shape as a function of 
the time r. Specifically, the frequency distribution 
of each coherence spectrum (Fourier transform on t) 
changes with the reversion time: higher frequency com¬ 
ponents decay faster. This distinctive behaviour, called 
eigen-selectivity [S], is the fingerprint of decoherence 
associated with an open quantum spin system. 


III. EXPERIMENTAL SETUP AND RESULTS 


In this section we present the experimental results on 
liquid crystal samples in the nematic phase, obtained 
with the pulse sequences shown in Figure [T| and de¬ 
scribed in Section ttn 

The experiments which use the reversion sequence 
MREV8 were carried out on samples of 5CB (4’-pentyl- 
4-biphenyl-carbonitrile) and PAAde (methyl deuterated 
p-azoxyanisole), using a home-built pulsed spectrometer, 
based on a Varian EM360 magnet (60MHz for protons). 
The electronic setup allows control of the pulse phases, 
homogeneity of the magnetic field, and sample temper¬ 
ature with stability of ±0.1°C. Experiments which use 
the reversion sequence MS were conducted on a Bruker 
minispec mq20 (20 MHz) on 8CB (4’-octyl-4-biphenyl- 
carbonitrile). 

The experiment starts by setting the time tp needed to 
prepare the initial state ps{tp) so that a subsequent free 
evolution drives the system to the so called “strong” {S- 
order) or “weak” W-order [7]. We test the experiment by 
preparing the initial state in 5CB with tp = 27/rs, since 
it has already been shown that after a transient t ~ 1ms 
the system reaches the S quasi-equilibrium state m, 
with only zero-quantum coherence (on the z-basis). Eig- 
ure[^ shows the experimental results on 5CB. The data 
is organized so that the /i-axis displays the spectra of 
the different coherence orders Uc (uc = 0, ±1, ±2,...). 
The spectrum of Uc order is symmetric to the one with 
—TT-c order, and the figure is arranged to show all the 
resolved coherence orders. The phase ip was varied in 
steps of 9.47 degrees in order to encode coherences with 
\nc\ < 19, however, spectra of coherences |nc| > 3 have 
negligible intensity and are not shown. The r-axis shows 
the spectra measured at increasing reversion times. In 
Figure (a) the reversion block is a single MREV se¬ 
quence where the total reversion time r = 12 ti is varied 
by increasing the delays ti (as shown in Figure[^a)). The 
largest contribution corresponds to zero-order coherence 


and, as expected, its spectrum is sine-shaped (because of 
inevitable truncation of a zero frequency signal). A much 
lower amplitude single-quantum spectrum, observable at 
short times is seen to attenuate for r ~ 400/iS. A similar 
behaviour is obtained in Figure]^ (b), where the reversion 
block is a series of concatenated MREV8 blocks with to¬ 
tal length r = nTc, where n is the number of blocks and 
the characteristic step-time is Tc = 12ri = 90.36p,s. This 
sequence yields better results since the zero-order coher¬ 
ence amplitude has a much slower attenuation with r, 
indicating a good reversion power, consistent with the 
ideal case calculated in the Appendix. 




'[ms] 


FIG. 2. Experimental spectra of the coherences measured in 
nematic 5CB. The time tp is set to prepare the strong or S- 
order {tp = 27ps). The /r-axis shows a combination of spectra 
obtained for different coherence orders, that are symmetric 
with respect to zero order. Each coherence spectrum is limited 
to ±25 kHz. The reversion block in (a) corresponds to a 
single MREV8 sequence where the total reversion time r = 
12ri is varied in steps of 120/rs by increasing the delays ri; 
in (b) we used an increasing number of short blocks with a 
characteristic step-time of = 90.36/ts. 
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A richer coherence spectrum is obtained by preparing 
the initial (W) state with tp = 69/iS. In this case the zero- 
, single- and double-quantum coherence spectra are well 
distinguished, as shown in Figure]^ The amplitude of the 
zero-quantum coherence has a different dependence with 
T under the two initial conditions. This behaviour is a 
consequence of the non-ideal performance of the reversion 
sequence, as shown in the Appendix [A} 



FIG. 3. Experimental spectra of the coherences measured in 
nematic 5CB. The time tp is set to prepare the weak or W- 
order {tp = 69/rs). The /r-axis shows a combination of spectra 
obtained for different coherence orders, that are symmetric 
with respect to zero order. Each coherence spectrum is lim¬ 
ited to ±25 kHz. The reversion block corresponds to a single 
MREV8 sequence where the total reversion time r = 12ri is 
varied in steps of 120/rs by increasing the delays ti. 

The evolution of the multi-spin single-quantum coher¬ 
ence is worth to be analyzed separately: Figure Qa) 
shows a stack plot of the single quantum spectra /r = I 
obtained under the initial condition W for different val¬ 
ues of r. Notice that, in terms of Eq. ( [4^ , the shape of 
these single quantum {vi = 1) spectra correspond to 

Ai Ai 

It is then expectable that they differ from the NMR 
spectrum, which corresponds to the single-spin, single¬ 
quantum coherence Ai = Ai = I. 

In order to display the frequency content as a func¬ 
tion of the reversion time, Figure |^b) shows several cuts 
of the amplitudes of the spectra at fixed frequency val¬ 
ues, normalized to the corresponding values at t = 0. It 
can be seen that the higher frequency components decay 
faster than the low frequency ones, and that the charac¬ 
teristic decay times are noticeably different, the one of the 
higher frequency being approximately one third of that of 
the lower frequency. Previous work on 5CB [5] and 8CB 


m gave a first indication that the attenuation the over¬ 
all amplitude of the reverted single-spin, single-quantum 
signals is a consequence of irreversible decoherence. The 
results presented confirm that and also show that the 
higher frequency components are selectively attenuated 
while the quasi-equilibrium is developing. 

The double quantum peaks have a much smaller am¬ 
plitude than the multi-spin single quantum ones, just as 
happens with the calculated ones (see Appendix A), and 
their S/N ratio is too small to analyze the selective nar¬ 
rowing. However it can be seen that their amplitude 
decays monotonously with r. It is worth to point out 
that this behaviour differs from the calculated one, since 
as shown in Figure |^c2), the calculated amplitude of 
coherence ±2 has the same trend as coherence 0. This 
difference indicates that the definite decay of this coher¬ 
ence observed in the experiment is due to a mechanism 
not accounted by the closed-system dynamics. 

The behaviour of a sample of PAA^^g under the same 
reversion sequence (concatenated blocks) is shown in Fig¬ 
ure]^ where (a) corresponds to a preparation time for the 
S condition and (b) to the W condition. The phase ip was 
varied in steps of 20.5 degrees. Again the W condition 
gives rise to a richest coherence content. The amplitude 
variation of the zero-order coherence in both plots (a) 
and (b) can be explained by the inherent non-ideality of 
the reversion sequence (see Appendix A). A closer anal¬ 
ysis of the selective attenuation of the multi-spin single¬ 
quantum spectra of PAA^^g is shown in Figure where 
(a) contains a superposition of all the spectra at different 
reversion times r, and (b) shows the r dependence of the 
amplitude of different frequency components of the spec¬ 
tra obtained from (a). The data plotted in (b) is normal¬ 
ized in order to spot on their characteristic decay times. 
The dashed lines show the exponential decay fitting to 
each frequency component decay curve. The character¬ 
istic times depend on the frequency, being shorter for 
the larger frequency components, similarly to the case of 
5CB. 

In order to show that this characteristic behaviour 
is independent of the particular reversion sequence and 
compound, we carried out a similar experiment using 
the reversion sequence known as Magic Sandwich (MS) 
instead of MREV on the nematic liquid crystal 8CB. 
Figure [^a) shows the multi-spin single-quantum spec¬ 
tra at different reversion times in 8CB, starting from 
the W condition {tp = 72/rs). In this reversion sequence 
r = 1.5 tm as in Figure [^c) pT] . 

The amplitude of the higher frequency components 
decrease with r faster than the central peak. Figure [^b) 
shows the r dependence of the amplitudes of the central 
and side peaks. The central peak amplitude decay was 
fitted with a linear function with slope = 3.84ms, 
while the characteristic decay time of the side peak 
(exponential fit) is = 1.24ms. 
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FIG. 4. (a) Front view of the spectra of coherence order /.i = 1 
obtained in Figure for nematic 5CB. (b) Time evolution 
of the (normalized) amplitude at different frequencies on the 
spectra of (a). The dashed lines are exponential fittings with 
different characteristic times as shown in the legend. 

IV. DISCUSSION AND CONCLUSIONS 

A wide variety of systems develop quasi-equilibrium 
states in an early time-scale, that is, in short times as 
compared with the lapse necessary to transfer a signif¬ 
icant amount of energy to the environment. Once the 
system attains a quasi-equilibrium state, the dynamics 
of its spectroscopic properties can be described by a den¬ 
sity operator that has the form of an expansion in quasi¬ 
invariants. There are many advantages in the use of a 
quasi-equilibrium description for theoretical as well as 
for practical analysis of the system dynamics. For ex¬ 
ample (!) representing the density operator in terms of 
a basis of quasi-invariants implies a very simple analyt¬ 
ical representation that can be written as a linear com¬ 
bination of a few quasi-invariant operators; {ii) the fact 
that quasi-invariants are immune to decoherence, turns 


FIG. 5. Experimental spectra of the coherences measured 
in nematic PAA^je at T = 115°C, using blocks of MREV8 se¬ 
quences of duration Tc = 88.36/rs, which define the r axis step- 
times. The time tp is fixed to set to prepare different dipolar 
quasi-equilibrium states, (a) Strong or 5-order {tp — 40/rs), 
(b) Weaker or W-order {tp = 84/is). The /r-axis shows a com¬ 
bination of spectra obtained for different coherence orders. In 
(a) each coherence spectrum is limited to ±5 kHz and in (b) 
to ±10 kHz. 


them into good candidates to be used as memory units 
in quantum computations or to preserve information of 
the dynamics, etc.;{iii) since quasi-invariants only evolve 
by relaxation processes, they are very useful as sensors of 
the molecular motions and other long time scale effects 
of the environment over the spin dynamics. Therefore, 
an adequate understanding of the processes that bring 
the system toward the quasi-equilibrium states becomes 
essential. 

In the literature on solid state NMR, the assumption is 
frequently made that the spin system can be considered 
as strictly closed in the early time-scale, and that the 
spin-lattice relaxation processes are the only channel by 
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xtMS] 


FIG. 6. Front view of the single quantum spectra of nematic 
PAAde prepared in the W condition (/r = 1 of Figure [^b)), 
corresponding to r < 1400/is.(b) Time evolution of the am¬ 
plitude of different frequency components of the spectra ob¬ 
tained from (a). The curves were normalized to compare their 
decay rates. The dashed lines are the exponential fittings to 
each frequency component decay. Label indicates the charac¬ 
teristic decay times r^. 


which the spin environment can act on the spin system 
HOI EH- In this vision, the quasi-equilibrium state is not 
an actual state, but merely an apparent representation 
of the dynamics produced by the effect of a phase super¬ 
position from the spin variables, when an observable is 
measured. Following this vision, the decay of the spin 
observables when reversion experiment are performed is 
attributed either to imperfections in the experimental 
settings of the pulse sequences, or to the dissipative in¬ 
teraction among all the spins in the sample. Ideally, after 
improving the experimental techniques, the only limit to 
reversion would be some kind of ‘internal equilibration’ or 
‘pseudothermalization’, a unitary process that would be 
the cause of the observed coherence decay Although this 


0.5 1.0 1.5 




FIG. 7. (a) Stack plot of the spectra of the single quantum 
coherence in 8GB for different reversion times r. (b) Time 
evolution of the (normalized) amplitude of the central peak 
(blue circles) and the side peaks (green triangles). The dashed 
lines are exponential fittings which show the marked difference 
in the characteristic decay times of the two frequency curves. 
The central peak amplitude decay was fitted with a linear 
function with slope Td = 3.84ms. 


line of thought might in principle give a phenomenologic 
explanation the observed signal attenuation in reversion 
experiments in solids, it cannot explain the phenomenon 
in liquid crystals. In the LC mesophases only the few 
spins which belong to each molecule are effectively inter¬ 
acting through their magnetic degrees of freedom, there¬ 
fore the interaction among a thermodynamic number of 
spins cannot be claimed as the mechanism which inhibits 
longer experimental reversion times. 

On the other hand, in the field of open quantum sys¬ 
tems, the adiabatic decoherence due to the quantum in¬ 
teraction of the observed spin system and the reservoir 
is shown to be responsible of the decay of the matrix el¬ 
ements of the density operator |M1 [H] . Particularly, in 
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references [S] and [5] the attenuation of single-spin NMR 
signals in reversion experiments on liquid crystals is ex¬ 
plained under the formalism of open quantum systems. 
Such attenuation is essentially irreversible, even in the 
assumption of a perfect-reversion experiment. The char¬ 
acteristic behaviour of decoherence in an open quantum 
system is eigen-selectivity, which produces the character¬ 
istic frequency compression on the spectra of the coher¬ 
ences. 

The importance of clarifying which class of spin system 
must be treated as an open quantum system in NMR ex¬ 
periments resides in the fact that spin systems are good 
candidates to be used as quantum registers, and in such 
case it is crucial to have longer decoherence times to al¬ 
low the execution of a larger amount of quantum gates. 
Therefore, it is very important to know which descrip¬ 
tion of the spin systems is more adequate, the closed 
or the open one, since it would lead to devising differ¬ 
ent encoding strategies. The results obtained in Sec¬ 
tion III contribute to this discussion. The experiment 
proposed in this work allowed monitoring the evolution 
of the multi-spin quantum coherences under a reversed 
spin dynamics. The single and double quantum coherece 


amplitudes attenuate monotonously under different re¬ 
version schemes in a way that definitely differs from the 
evolution of a closed system under non-ideal experimen¬ 
tal conditions. Besides, the frequency content of the 
multi-spin single-quantum coherence of the three stud¬ 
ied compounds changes with the reversion time, showing 
that the higher frequency components attenuate faster 
than the lower ones. We interpret this results as an 
evidence of the occurrence of the eigen-selectivity effect 
during the transient from an initial state to the quasi¬ 
equilibrium state. This poses decoherence as the mech¬ 
anism which drives a system to a QE state. Thus, we 
can conclude that the quasi-equilibrium states are actual 
states that the open quantum systems may attain during 
their transit to thermal equilibrium. 
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Appendix A: Analysis of non-idealities in the 
spin-dynamics reversion sequence. 

This appendix is dedicated to analyze the effects of 
non-idealities in the dynamics of the spins under the 
reversion sequence MREV8. In particular, we address 
the contribution of the non-secular terms in the dipolar 
Hamiltonian, which generally emerge when a period of 
unitary evolution under the dipolar Hamiltonian is en¬ 
closed between two pulses, as happens in the reversion 
sequences. 

To mitigate the effects of the non-secular part over 
the dynamics, the free evolution time is usually set as 
small as possible, but not shorter than the irradiation 
time. However, this strategy does not apply when long 
reversion periods are needed. 

We performed a numerical calculation of the signals 
produced by the experiment described in Section |llj on 


an 8-spin system having the geometry and dipolar cou¬ 
plings of the ^H nuclei in a PAA^g molecule [8]. The dy¬ 
namics was calculated considering the spins as a closed 
system. In order to expose the effects of the non-secular 
dipolar terms, we calculated the outcome signals when 
the reversion sequence is a series of MREV8 blocks, we 
explored two different settings of the characteristic length 
Tc of the individual blocks and their expected effect over 
different values of the preparation time tp. 

The graphs on left column of Figurej^show the spectra 
corresponding to the condition ri = 5/iS and the left col¬ 
umn to Ti = 20/rs. The first row corresponds to tp = 0, 
or equivalently to reverting the FID (single-spin single¬ 
quantum signal); rows (b) and (c) correspond to prepa¬ 
ration of the 5-order, and W-order, respectively, and the 
preparation time in row (d) is an intermediate value be¬ 
tween the former. The total reversion time in the r axis 
is determined by the number of blocks. 

Reversion is optimal when the characteristic time is 
small (ti = 5/rs): the amplitude of each coherence 
is practically independent of r throughout the range 
scanned (I ms) in this simulation, as seen in Figure|^al), 
(bl), (cl) and (dl). On the contrary, a larger setting of 
Ti introduces variations on the amplitudes of the differ¬ 
ent spectra with the reversion time r, as shown in (a2), 
(b2), (c2) and (d2). Particularly, the amplitude of the 
zero-quantum coherence peaks in (b2) and (d2), calcu¬ 
lated in the closed system, show a slow attenuation simi¬ 
lar to the one observed in Figure [^a) and (b), as well as 
in Figure [^a). When the initial state is prepared with 
tp = 97.5^s, as in Figure [^c2), the amplitude of zero- 
quantum coherence decreases and grows similarly to the 
experimental spectra of Figure and that of Figure 
The analysis indicates that this feature on the evolution 
of the zero-quantum coherence is just introduced by the 
non-reverted terms of pulse sequence and not due to the 
open character of the evolution. 
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(a1) (a2) 



FIG. 8. Calculation of the proposed experiment for different 
preparation time, tp, of the state of a PAA^jg molecule (dipole 
couplings from Ref. [33]). al,a2: FID {tp = Ofj,s), al: ti = 
5fis, a2: t\ = 5/iS. bl,b2: 5-order {tp = 47.5/is), bl: ri = 
5/rs, b2: n = 20/is. cl,c2: W-order {tp = 92.5^s), cl: n = 
5/rs, c2: ri = 20fj,s. dl,d2: Mixed-state {tp = 195/rs), dl: 
Ti = 5/iS, d2: ri = 20/iS. As a guide for visualizing the 
evolution over the time r, solid lines are traced joining the 
maximum peak of the spectra of different coherences {fj, — 0 
in red, /r = 1 in blue, ^ = 2 in green, jj, — 3 in magenta). 





















